0.4 




0.75 0.80 0.85 0.90 

n 



0.6 I— 



+ 



0.4 - 



0.2 - 



0.0 - 




MMP model 
Hubbard model 
harmonic 
anharmonic 
(Y,Pr)Ba2Cu30, 



-0.2 



25 



75 
T. [K] 



125 



THE ISOTOPE EFFECT 
IN d-WAVE SUPERCONDUCTORS 

H.-B. Schiittler and C.-H. Pao 
Center for Simulational Physics 
Department of Physics and Astronomy 
University of Georgia 
Athens, Georgia 30602 

Abstract 

Based on recently proposed anti-ferromagnetic spin fluctuation exchange models 
for c/j;2_j^2 -superconductors, we show that coupling to harmonic phonons cannot 
account for the observed isotope effect in the cuprate high- Tc materials, whereas 
coupling to strongly anharmonic multiple-well lattice tunneling modes can. Our 
results thus point towards a strongly enhanced effective electron-phonon coupling 
and a possible break-down of Migdal-Eliashberg theory in the cuprates. 
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A growing, but not uncontroverted body of experimental evidence currently 
suggests that Y Ba^Cu'^O^ and, possibly, other cuprates are c/j,2_j^2 -super- 
conductors."^'^ On the theoretical side, (/-wave superconductivity would probably 
rule out conventional phonon-mediated pairing,^ while supporting an electronic 
pairing mechanism, most prominently anti-ferromagnetic (AF) spin fluctuation 
exchange models. ^'"^'^ Yet, except for certain "optimal" doping concentrations, 
many cuprates, including YBa2Cus07, exhibit a quite noticeable isotope effect,^ 
the classical hallmark of the superconducting electrons' coupling to the lattice 
vibrational degrees of freedom. 

Starting from recent diagrammatic AF spin fluctuation models, ^'"^'^ we will 
show here that the observed order of magnitude of the isotope exponent^ a ~ 1 does 
in fact imply an extraordinarily large electron-phonon (EP) coupling parameter 

1. For coupling to harmonic phonon systems, such coupling strengths would 
be well in excess of upper limits imposed by structural stability arguments and, 
independently, by the observed Tc and normal-state transport scattering rates. 
However, the required A's are consistent with a coupling to local, large-amplitude 
lattice tunneling excitations, arising in a strongly anharmonic multiple-well lattice 
potential. Barring experimental problems in the isotope measurements,^ our re- 
sults thus point towards a very strong enhancement of the effective EP coupling, 
strong enough, in fact, to cause a break-down of the very approximation upon 
which the conventional Migdal diagrammatic theory is based. As such, our results 
point towards the possibility that lattice vibrational degrees of freedom play a 
central role in the low-energy electronic properties of the cuprates. 

We start from the linearized Migdal- Eliashberg equations^ for the self-energy 
E and pair wavefunction $ at wavevector k and Matsubara frequency iv = {2m -\- 
l)7rzT (in energy units) at the transition temperature T = Tc, 

E(A;', iv') = -j^^y-{k' -k, iv' - iv)G{k, iv) 

J. (1) 

$(A;',z//') = ~ ~ iT^)\G{k,iu)\^^k,iu). 
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Here, k = (kx^ ky) is summed over an N = LxL grid covering the Brillouin zone of 
a two-dimensional (2D) square lattice. The single-electron Green's function G is 
obtained self-consistently from G(A;, zi^) = [iu — ek — 'S{k,iu)]~^ , assuming a single, 
2D electron band €k = —2ti(cos(kx)+cos(ky))—4:t2Cos(kx)cos(ky) — fi with chemical 
potential ji and 1st and 2nd neighbor hopping ti (= 250meV)^" and t2 , respectively. 
For singlet d-wave pairing, the effective electron-electron interaction potentials in 
Eqs.(l) are V±{q,iuj) = ib(/^Xs(g, zcu) -|- Vp{q,iuj). This includes a spin-fluctuation 
term with coupling constant and dynamical spin susceptibility Xsi<l, ito),'^''^'^ and 
a phonon term of the general form Vp{q,iuj) = —Upf{q)Q'^/{uj'^ -\- fig), assuming 
a single phonon branch with an energy dispersion Og and a form factor f{q)^ 
normalized so that Up = —N~-^ Sg^l^'O)- ratio A = Up/B provides then 

a rough estimate for the dimensionless Eliashberg parameter A,^ averaged over 
the electronic bandwidth B = 8ti. We emphasize that the isotopic mass- (M-) 
dependence of Vp enters only via Og, but not via Up.^ 

Global structural stability of the undoped, ^-flUed-band cuprate parent com- 
pounds requires the phonon-mediated on-site attraction Up to be less than the 
Coulombic on-site repulsion, of order of the on-site Hubbard ?7-parameter.^ If 
Up > U ^ the conduction electrons would form local pairs (bi-polarons), rather than 
local magnetic moments. At |-fllling, such local pairs would undergo ordering in 
a charge superlattice,^ rather than forming the AF spin superlattice observed"^ 
in the undoped cuprates. This transition from local moment AF spin density or- 
der to local pair charge density order has been established in recent studies of 
Holstein-Hubbard and related EP models.^ 

Local structural stability in the doped systems, i.e. stability of the dopant in- 
duced carriers against polaronic self-localization,^ requires an even more stringent 
upper hmit Up < U^p°"^ < U. If Up > U^p°"\ pol aron formation would destroy (or, 
at the very least, strongly renormalize and broaden) the delocalized quasi-particle 
states whose existence is essential for the Migdal- Eliashberg approach. Holstein- 
Hubbard estimates give 

^{loc) _ 2 - 3 X ti for a Hubbard ~ 8 - 12 x t 1 in the 
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cuprates/ or, equivalently, A('°") = ^/^"'V^ ~ 0.25-0.4 and A^^'"^) = ~ 1 - 1.5. 

Our results for the dj.2_y2 superconducting transition temperature Tc and its 
isotope exponent a = —d\og{Tc) j d\og{M) are illustrated by Figs. 1 and 2 for the 
case of an Einstein phonon model with Og = Oq = const, local coupling f{q) = 1, 
and assumed isotopic mass dependence Oq M~^/'^ . Such a model could roughly 
represent, for example, the local coupling(s) to high-energy optical modes which 
are dominated primarily by the lightest atomic species, i.e. oxygen, in the cuprates 
with Oo ~ lOOmeV.^ The Xs in the MMP model^, used in Figs. 1 and 2, is given 
by 

Imx.(<Z,^ + «0+) = Im[xQ/(l + e'k-QI' -«^M)]%c - \io\) (2) 

with Q = (7r,7r) for qx^Qy > 0. Except for values explicitly stated in Figs. 1 and 
2, we have used Up = leV, Oq = 50meV and the band and T-independent spin 
fluctuation parameters from Table II of Ref. [5a], hereafter referred to as MP-II, 
for a hole doping concentration x = 1 — n = 25%. Eqs. (1) are solved accurately, 
without further approximations, by fast Fourier transform techniques."^" 

Fig. 1 shows Tc and a as functions of the electron concentration n. To model 
roughly the doping dependence of x^, iOs{x) and S,{x) are varied with x = 1 — n. 
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like cUsOCi^" OCX, and attain their MP-II- values at a; = 25%. Fig. 2 shows a 
Tc for the same model at flxed a; = 25%, with Tc varied by increasing Up from to 
2eV. Fig. 2 also shows results for coupling to anharmonic phonon and a Hubbard- 
based spin fluctuation model, discussed below. The flgures clearly demonstrate the 
very general result that coupling to harmonic phonons will substantially suppress 
the d-wave Tc without, however, causing a large enough isotope effect to account 
for the order of magnitude of the experimental results, illustrated here by the 
{Y,Pr)Ba2Cu307 data.^" 

Note here that the MP-II spin fluctuation parameters, with gg = 640meV, were 
carefully tuned^" to reproduce the observed NMR data, resistivity and Tc = 90 °K 
in Y Ba2Cus07 ^ excluding the effect of phonons. If we now include phonons, with 
a strength of only Up = leY (i.e. A = 0.5), say, we need to roughly double g^ (see 



Fig. 1) in order to maintain a Tc = 90 °K. Thus, we will also double, roughly, the 
normal-state model resistivity from spin fluctuation scattering (not even counting 
phonon scattering !), which already seriously compromises the agreement with the 
experimental resistivity data. So, in the MMP model,^ Eq. (2), the observed 
maximum Tc = 90 'K, combined with the resistivity data for Y Ba2Cus07 ^ will 
again limit Up to Up^leV. 

We have generated overwhelming numerical evidence that the foregoing results 
are generic and robust, (i) by systematically varying, over wide ranges, all spin 
fluctuation, electron-phonon and band parameters (ii) by using different types of 
form factors f{q)^ different types of and/or multiple g-dependent phonon branches 
(including e.g. 3D acoustic phonons), (iii) by modifying the AF spin fluctuation 
model (introducing e.g. spin gap effects into the MMP model or using other pro- 
posed types of spectral models^'"^"^'"^"^), and (iv) by adding impurity scattering. "^"^'"^^ 
While such modiflcations of the model may substantially affect Tc, they do not 
increase the overall magnitude of a signiflcantly. We thus conclude that the ob- 
served magnitudes of Tc and a in the cuprates cannot possibly be explained by 
a conventional Migdal-Eliashberg treatment of the conduction electrons' coupling 
to AF spin fluctuations and harmonic phonons. 

To understand why harmonic phonon exchange Vp causes a large Tc suppres- 
sion, but only a ^ma// isotope effect, consider flrst its effect on the self-energy E in a 
simple McMillan-type analysis, ^'"^ generalized to the case of d-wave pairing. In the 
real-frequency domain, iv^v -\- iO~^ , Vp contributes, via V-, to the quasi-particle 
mass enhancement Z = 1 — c?i/ReE|j^ = 1 -\- Zg -\- Zp and to the quasi-particle 
damping ImE = ImEs -|- ImE^, with "s" and "p" denoting the spin fluctuation and 
phonon contribution, respectively. Now, Zp = X^Up/ B is independent of isotopic 
mass and, in McMillan theory, enters into Tc roughly as^'"^ Tc ~ {u) exp( — Z / Xd) 
with a bare d-wave pairing strength^'"^ Xd and an appropriate average energy {u) 
of the V^-spectrum.^'"^ Thus, Zp causes a large isotope-independent reduction in 
Tc, by roughly a factor exp( — A/Aj/), which is quite consistent with our numerical 

- 5 - 



results for the ?7p-dependence of Tc- On the other hand, the phonon contribution 
to the quasi-particle damping,^ |ImEp| ~ AmaXg(Og), is comparable to the phonon 
energy scale Og and thus much smaller than the typically electronic magnitude 
of E^. Therefore, ImE^ causes an additional Tc-suppression, by pair-breaking, 
which does depend on isotopic mass, but is much smaller than the "intrinsic" 
Tc-suppression, due to E^, thus causing only a small a.^^ Treating the phonons 
solely in terms of an Abrikosov-Gorkov pair-breaking theory-^"^ (without inclusion 
of the mass enhancement effect !) is seriously flawed, since it severely underes- 
timates the Tc-suppression caused by the phonons, for given magnitude of the 
isotope exponent. In this regard, the effect of phonons is quite different from that 
of impurities. "^"^'"^^ 

In the case of a g-dependent V^(g,zc<;), Vp may also contribute directly, via 
V+, to the (/-wave pairing potential. In McMillan theory,^ '"^ the dominant effect is 
again an isotope-independent phonon contribution Xd^p to the bare d-wave pairing 
strength Xd = Xd^s + Xd^p entering into Tc ~ (u) exp( — Z/Xd). The isotope-dependent 
contribution {io)p to (cu) = {io)s + (^)p is again of order of the phonon energy scale, 
i.e. small compared to the spin fluctuation contribution (cu)^, thus, again, giving 
only a small contribution to a. A substantially larger a and smaller Tc-suppression 
may be obtained in the case of phonon-mediated d-wave pairing, i.e. if Vp(q^iiLi) 
is attractive in the c/j,2_j^2 -pairing channel, with large enough Xd,p>0. Because 
of the expected cancellation effects in the presence of multiple phonon branches, 
this scenario seems unlikely, but should nevertheless be explored by flrst-principles 
Ad,p-calculations. 

Phonon renormalizations of the spin fluctuation exchange potential g^Xs do 
also not provide a viable mechanism for obtaining a signiflcantly larger as 
illustrated by the Hubbard model results in Fig. 2. Here, both E and g^Xs, as well 
as a less important charge fluctuation term g'^Xc in V±, Eq. (1), were calculated 
self-consistently, from G{k,iv), in the fluctuation exchange approximation to the 
2D Hubbard model, ^'"^^ using numerical renormalization group techniques. The 
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phonon term Vp in Eq. (1) thus explicitly modifies Xs Xc, via E. The result, 
|a| ^ 1, is explained by noting that, again, the dominant effect of Vp on Xs ^-nd Xc 
comes from the isotope-independent mass renormalization Zp. 

Fig. 2 also shows results obtained for coupling to an anharmonic phonon 
system, consisting of independent, local anharmonic oscillators with one lattice 
displacement degree of freedom uj per lattice site j. Each uj is subject to a 
local Hamiltonian hj = —{Ti'^/2M)d'^/du^- + w{uj) with a double-well potential 
w(uj) = AAsi^iuj / dY — {uj /dY]. Here, d is the distance between the two local w- 
minima. The chosen tunneling barrier height = 26.89meV and re-scaled atomic 
mass M = Md^ITi = 0.3555meV"^ give a double-well tunneling splitting = 
El — Eo = 7.25meV and a quasi-harmonic (single- well "phonon") excitation energy 
0/j = E2 — Eo = 35.5meV where Ei denotes the ^-th excited state energy of hj. The 
g-independent Vp was calculated as"^^ Vp{iijj) = —C^ j'l^ dr (uj{T)uj{0)) exp{iijJT), 
with the local correlation function {uj(T)uj(0)) obtained from numerical solutions 
for the eigenstates of hj. Tc in Fig. 2 was varied by raising the re-scaled cou- 
pling C = Cd from to^'^'^^'* of 200meV and SOOmeV for the MMP^ [Eq.(2)] 
and Hubbard"^ model, respectively. In contrast to harmonic phonon models, 
\ = —Vp{0)/Br^C'^/{Bnt), at r = o,^"'-^^" is strongly enhanced here and strongly 
isotope-dependent, since the tunneling splitting varies exponentially with the 
isotopic mass M. For given Tc, this model thus leads to a much larger a and, 
for given a, requires a much smaller bare coupling C than the harmonic phonon 
models discussed above. These larger values of a can be achieved without sub- 
stantially increasing the model's DC resistivity, since the phonon contribution to 
the quasi-particle damping, ImE, is, again, small compared to the spin fluctua- 
tion contribution. When studied as a function of doping concentration x, with 
1^ and/or iOa varied with x as in Fig. 2, the isotope exponent in the anharmonic 
model typically increases in magnitude with decreasing Tc as one moves away from 
the optimal doping. This behavior appears to be generic to the model and is in 
qualitative agreement with the data in the cuprates. 
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Substantial experimental evidence for a strongly anharmonic multiple-well 
lattice dynamics in the cuprates exists.^ Polaron formation is a possible mecha- 
nism for generating anharmonic tunneling modes of sufficient effective coupling 
strength to the conduction electron system within physically reasonable parame- 
ter limits.^ The foregoing treatment extends the Migdal-Eliashberg approach to 
model the two essential features of such a polaronic system, namely, multiple-well 
lattice anharmonicity and isotope-dependent electron mass renormalization.^'-^^" 
However, as polaronic instabilities may well cause a break-down of the Migdal 
approximation, the present treatment should be re-examined in a more general 
strong-coupling framework.^ 

In summary, based on a McMillan analysis and large body of numerical evi- 
dence, we conclude that harmonic phonon exchange causes a large Tc-suppression, 
but only a small isotope exponent a, in d-wave pairing instabilities driven by 
AF spin-fluctuation exchange. Physical constraints on the EP coupling strengths 
limit |a| to theoretical values below 0.1, which is an order of magnitude smaller 
than observed values in non-optimally doped, reduced- Tc cuprate materials. By 
contrast, acceptable values of a at reasonably large T^s can be obtained, with 
physically reasonable coupling constants, by exchange of strongly anharmonic lat- 
tice tunneling excitations. Such anharmonic lattice modes lead to a very strongly 
enhanced, isotope-dependent effective A-parameter. As such, our results, com- 
bined with the experimental isotope data,^ can be regarded as evidence for a very 
strongly enhanced effective EP coupling in the cuprates which may well cause 
a break-down of the Migdal approximation. The further development of strong- 
coupling approaches,^ beyond Migdal-Eliashberg theory, will be crucial for gaining 
a deeper understanding of the role of phonons in the low-energy physics of the 
cuprate high- Tc systems. Experimentally, it would be of great interest to study 
the isotope-dependence of other physical properties, for example the speciflc heat, 
in order to test whether the cuprates indeed exhibit a strong isotope-dependence 
of their electronic quasi-particle mass, as predicted by polaronic and anharmonic 
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lattice models. 
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FIGURE CAPTIONS 



Fig.l. dj.2_y2 transition temperature Tc and isotope exponent a vs. electron 
concentration n for coupling to Einstein phonons in the MMP model, Eq.(2), 
with ujs oc oc (1 — n). Also shown are the Tc and oxygen isotope data for 
(Y^ Pr)Ba2Cus07 ^ from Ref. [6a], with n = 0.75 assumed at the maximum Tc- 

Fig. 2. Isotope exponent a vs. dj.2_y2 transition temperature Tc in the Hub- 
bard model, Ref. [4a, b], with U jti = 6, ti = 250meV, t2 = 0, and n = 0.86, and 
in the MMP model, Eq.(2), with MP-II parameters, for coupling to harmonic 
Einstein phonons, and for coupling to anharmonic tunneling excitations, with 
EP parameters given in text. Also shown are the oxygen-isotope-vs. -Tc data for 
{Y,Pr)Ba2Cu307, from Ref. [6a]. 
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